Bulk Fermi surface coexistence with Dirac surface state in Bi 2 Se 3 : a comparison of 
photoemission and Shubnikov-de Haas measurements 
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Shubnikov de Haas (SdH) oscillations and Angle Resolved PhotoEmission Spectroscopy (ARPES) 
are used to probe the Fermi surface of single crystals of Bi2Se3. We find that SdH and ARPES 
probes quantitatively agree on measurements of the effective mass and bulk band dispersion. In 
high carrier density samples, the two probes also agree in the exact position of the Fermi level Ef, 
but for lower carrier density samples discrepancies emerge in the position of Ef- In particular, SdH 
reveals a bulk three-dimensional Fermi surface for samples with carrier densities as low as 10 17 cm~ 3 . 
We suggest a simple mechanism to explain these differences and discuss consequences for existing 
and future transport studies of topological insulators. 
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Recently, a new state of matter, known as a topolog- 
ical insulator, has been predicted to exist in a number 
of materials: Bi^SR,, Bi 2 Se 3 , Bi 2 Te 3 and Sb 2 Te 3 [H2|. 
This state of matter is characterized by a full band gap 
in the bulk of the material, but with a gapless, dissipa- 
tionless surface state. The surface state is comprised of 
counter-propagating spin states, which create a disper- 
sion of a single, massless Dirac cone that is protected by 
time-reversal symmetry. The experimental realization of 
this state could mean significant advances in spintronic 
devices, quantum computation and much more besides. 

As a result there has been great excitement in the last 
year after the discoveries of various ARPES experiments 
[3|— 15|] and more recently from scanning-tunneling mea- 
surements 041] that such a state appears to exist in 
nature. Amidst this flurry of recent results, it is easy 
to forget that these same materials have been the sub- 
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FIG. 1: (a)Temperature dependence of two typical samples of 
Bi2Se3 with carrier densities differing by two orders of mag- 
nitude, (b) Shows the carrier density n e , (c) the resistivity po 
at T=2K and (d) the mobility for samples of different thick- 
nesses. Each sample was a cleave from a parent sample, so 
that the surface area of each sample was kept constant. 



ject of careful and thorough research for much of the 
20th century. However, common to all the unambiguous 
measurements of the Dirac cone is the use of surface- 
sensitive probes. Only recently have transport measure- 
ments e merg ed specifically investigating the surface state 
(Refs [S|-lri(), all of which note the dominance of the 
bulk conductivity. It is thus of great interest to perform 
a coordinated study of these materials using both bulk 
transport experiments and surface sensitive ARPES ex- 
periments. Here we report results of these investigations. 
The transport experiments reveal quantum oscillations 
that indicate a bulk band structure and Fermi surface 
volume that monotonically change with doping. For car- 
rier densities in the range ~10 19 cm -3 , the transport ex- 
tracted band structure is in quantitative agreement with 
the bulk band structure determined by ARPES which 
also observes the Dirac dispersion of the surface state. 
The quantitative agreement between ARPES and SdH 
provides additional support for the existence of novel 
band structure in these materials. For lower carrier den- 
sity samples down to 10 17 cm -3 we observe SdH oscil- 
lations which unambiguously pin the Fermi level in the 
bulk conduction band, with a high level of consistency 
across all samples measured from the same batch. While 
ARPES places Ep near the SdH level for some samples, 
there are others from the same batch whose Ep is found 
to reside into the bulk gap. We discuss possible expla- 
nations for these discrepancies and the implications for 
transport studies of surface Dirac Fermions in samples 
near a metal-insulator transition. 

The material Bi2Se3 can be grown without the intro- 
duction of foreign dopants as either n or p type 13, 13] 
though is more commonly found as the former because 
the dominant defects tend to be Se vacancies. Quantum 
oscillatory phenomena, which provides evidence of bulk 



metallic behavior has been reported by Kohler et al. [12 1 
on low carrier density samples and more recently by Kul- 
bachinskii ,14] on high carrier density samples. Below a 
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FIG. 2: (a) ARPES band dispersion on samples of Bi2Se3 
with carrier density 2.3xl0 19 cm" 3 (batch S4). (b) Due to 
the quantization of the energy spectrum into Landau levels 
(LLs), oscillations appear in the magnetoresistance known as 
SdH oscillations. The SdH oscillations here are for a sample 
taken from the same batch as in (a) at 6 = 0, corresponding 
to a oscillatory frequency of F= 155T, consistent with Ef ~ 
160meV. ARPES and SdH are in good agreement for these 
high carrier density samples. 



carrier density of 7x 10 18 cm -3 , the band structure is well 
approximated by a single parabolic band, making the in- 
terpretation of transport measurements transparent [lij . 
Two n-type samples with carrier densities differing by 
two orders of magnitude are shown in Figure [1] For the 
low carrier density samples an upturn in the resistivity 
is seen, which levels off at sufficiently low temperature. 
This behavior has been attributed to the presence of an 
impurity band whose thermally activated conductivity is 
comparable to the band conductivity until carriers freeze 
out at around 30K [HI EH- This behavior is not ap- 
parent in the higher carrier density materials, where the 
band conductivity always dominates. Even though we 
have reduced the carrier density by 2 orders of magni- 
tude, the resistivity increases by one, suggesting that the 
mobility has increased in the low carrier density samples, 
consistent with previous measurements [13] • The low car- 
rier density samples are around an order of magnitude 
smaller than those of reported topological insulators in- 
cluding Sn doped Bi 2 Te 3 (n e ~ 1.7 x 10 18 cm- 3 ) Q or Ca 
doped Bi2Se3(n e ~5x 10 18 cm _3 )[5j and as a result may 
be better candidates in which to observe the transport 
properties dominated by the topological surface state. 

Single crystals of Bi 2 Se 3 have been grown by slow cool- 
ing a binary melt. Elemental Bi and Se were mixed in 
alumina crucibles in a molar ratio of 35:65 for batch SI 
(n e = 5x 10 17 ), 34:66 for batch S2 (n e = 3 x 10 17 ), 34:66 
for batch S3 (n e = 2.3 x 10 17 ), and 40:60 for batch S4 
(n e = 2.3 x 10 19 ). The mixtures were sealed in quartz 
ampules and raised to 750 °C and cooled slowly to 550 
°C, then annealed for an extended period. Crystals can 
be cleaved very easily perpendicular to the (0 1) axis. 
Measurements of the resistivity and Hall effect were mea- 
sured in a 14T PPMS using a standard 4-probe contact 
configuration and Hall measurements were performed us- 
ing a 6-probe configuration. For the latter, only data 



which was linear in the low field limit was used to avoid 
mixing with longitudinal components. In addition to this 
precaution, signal from positive and negative field sweeps 
was subtracted to extract the odd (Hall) components of 
the signal, after which the carrier density is extracted in 
the usual way. ARPES measurements were performed at 
beam line 10.0.1 of the Advanced Light Source (ALS) at 
Lawrence Berkeley National Laboratory. Measurement 
pressure was kept <3x 10 11 Torr, and data were recorded 
by Scienta R4000 analyzers at 15K sample temperature. 
The total convolved energy and angle resolutions were 
16meV and 0.2° (i.e. KO.OO?^ 1 ) or <0.012(ll- 1 ) for 
photoelectrons generated by 48eV photons) , at which en- 
ergy the cross-section for both surface state and bulk 
bands is strong. 

In Figure [2] we show complimentary ARPES and SdH 
data on samples from the same batch, with carrier den- 
sity determined by the Hall effect of n e — 2.3 x 10 19 cm~ 3 . 
The SdH reveals an anisotropic pocket of frequency 155 
T, corresponding to a filling of around 160meV (the band 
structure is not parabolic at this filling and so we assume 
similar band structure parameters as Kohler et aZ. [T3| 
characterizing similar carrier density samples of Bi2Ses). 
ARPES results on samples from the same batch, show 
the Fermi level 150meV above the bottom of the conduc- 
tion band in good quantitative agreement. Similarly, the 
effective mass (see below) extracted from SdH is in good 
quantitative agreement with that measured by ARPES. 

In Figure [3] we illustrate angle dependent SdH data 
(a) taken at 1.8K, on a sample from batch SI with a 
lower carrier density of 10 17 cm -3 . The SdH signal re- 
veals a pocket that is approximately an ellipsoid elon- 
gated about the c3 axis, consistent with measurements 
by Kohler et al. from the 1970's[l2j on samples with 
similar carrier densities. For a two-dimensional pocket 
expected from the surface state, quantum oscillations 
should vary as l/cos#, where 9 is the angle between the 
c3 axis and the field direction, so the present observations 
must originate from a 3D Fermi surface existing in the 
bulk. It has been shown by Kohler et al. and more re- 
cently by Kulbachinskii et al. that the conduction band 
structure for these low carrier densities is approximately 
parabolic L2J, [l4( , and so the band filling can be estimated 
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in Fourier space. We estimate the Fermi energy to be 
18meV above the bottom of the conduction band. 

In Figure |3] (b) we show the derivative of the lon- 
gitudinal magnetoresistance of a sample from batch S3 
and a fit of the entire data set using the usual Lifshitz- 
Kosevich formalism, to extract the effective mass and 
Dingle temperature To, with fit shown in (c). Fitting 
the entire data set, which is often more accurate than 
tracing the amplitude of the Fourier transform, our fit 
yields to* = 0.15m e and To — 3.5K, for this frequency 
(F=14T). Similar data for samples from batch S4 give 
to* = 0.125m e , T D = AK and F=155T. The mean free 
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FIG. 3: (a) Magnetotransport for samples from SI. As the 
angle is swept the frequency of the oscillation varies according 
to the topology of the Fermi surface. For a two-dimensional 
pocket the expected dependence is 1/cosS (shown in green in 
the inset). The observed angle dependence is clear evidence 
for a closed ellipsoidal Fermi surface pocket, similar to that 
observed by Kohler et al. • Similar SdH data was gathered 
on batch S2 and S3 on a number of samples. Samples from 
batch S3 showing the temperature dependence of the deriva- 
tive of the SdH signal in (b) and a fit to the data shown in 
(c) from which the effective mass, Dingle temperature and 
oscillatory frequency can be extracted. 



path is calculated using the orbitally averaged velocity 
and scattering time extracted from Trj yields Isa ~ 60nm 
and ls4 ~ 220nm. This data is wholly consistent with the 
very complete SdH studies of Kohler et al. [H, EH and 
more recently by Kulbachinskii et al. 



14j . In addition, 



the data was reproduced with high consistency on a num- 
ber of samples from the same batch, and even on samples 
from different batches with similar growth parameters. 

ARPES data on samples from the same batches as 
those shown in Figure [31 determining the effective mass 
as m* — 0.13 in very good quantitative agreement with 
SdH. However, the exact placement of the Fermi level in 
the band structure reveals some disagreement. In Fig- 
ure 2] we illustrate photoemission data for two separate 
samples from batch SI. In (a) the Fermi level is near the 
bottom of the conduction band in agreement with SdH, 
while in (b) it is in the gap (about 60meV below the con- 
duction band), crossing the Dirac cone with apparently 
no bulk contribution. A number of samples from similar 
batches, such as batch S2 and S3, also have shown similar 
variation in Ep crossing the gap on some samples. While 
Ep determined by photoemission appears to show some 
variation, it is important to note that the measured Ep 
from ARPES is either near or below the SdH Ep. 

Despite the good agreement of SdH and ARPES on 
high carrier density samples, and the agreement of the 
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FIG. 4: ARPES data on samples of Bi 2 Se 3 from batch SI. 
The horizontal lines show the crossing of the Fermi level 
(E — Ep = 0) and the Dirac crossing, (a) Band structure mea- 
sured by ARPES results on samples from batch SI showed the 
Fermi level near the SdH level of ~ 15meV from the bottom 
of the conduction band, (b) Measurement on another sample 
from batch SI showed the Fermi level in gap. Some other sam- 
ples from S2 and S3 also showed the Fermi level in the bulk 
gap. The variation might be due to the lower carrier density 
of these samples, and the surface band structure is more sus- 
ceptible to small amounts of surface contamination, c) Band 
structure of a sample also from batch SI which was cleaved in 
atmosphere and exposed for 10s, showing significant n-type 
doping with large bulk conduction band pocket. 



effective mass and other band parameters on the low car- 
rier density samples, the discrepancy in the position of 
the Fermi level requires explanation. Such differences can 
occur for a number of reasons, for example due to sample 
variation within a batch, or perhaps due to variation in 
the exposure of cleaved surfaces before a photoemission 
measurement. However, it should be noted that the SdH 
frequency does not appear to vary significantly within a 
batch for up to 20 samples measured in the present study 
and so the former seems an unlikely scenario. Another 
reason for the discrepancy may be that atmospheric ex- 
posure of transport samples has contaminated them with 
an n-type dopant causing them to appear bulk n-type. 
Figure 2] (c) illustrates photoemission data for a sam- 
ple cleaved in air. The Dirac cone of the surface state 
remains robust and the bulk conduction band appears 
partially occupied. Such doping may lead to a 3D Fermi 
surface pocket appearing in SdH oscillations if the con- 
tamination is deep enough and allows for sufficiently long 
mean free paths. 

To investigate this possibility further we measure the 
thickness dependence of the transport by systematically 
thinning a single sample. Cleaving was achieved with 
tape, keeping the surface area of the resulting samples rel- 
atively constant and allowing direct comparison of data 
sets of each cleave. Though the samples are vulnerable 
to deformations between cleaves, only data from mirror- 
like flat samples is presented. In most cases these sam- 
ples still exhibited quantum oscillatory phenomena, con- 
firming the high quality of the cleaved samples. Figure 
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FIG. 5: (a) A schematic representation of the band structure 
seen by ARPES (solid) red horizontal line denoting the Fermi 
level as seen by SdH. (b) We infer band bending of about 
60 meV at the surface from a comparison of ARPES and 
quantum oscillations. 



[T] (b-d) shows a summary of the low-temperature car- 
rier density, resistivity and mobility. Within our error 
bars, each quantity seems to vary weakly down to 3/zm 
in thickness. Although the residual resistivity and carrier 
density varies slightly (possibly from disorder related to 
slight sample deformation, despite the precautions men- 
tioned above), the mobility remains almost constant as 
a function of thickness at /i ~ lm 2 /Vs. In summary, 
the transport is insensitive to the thickness, suggesting 
that the SdH oscillations are not a consequence of atmo- 
spheric contamination, but originate from the intrinsic 
band conductivity of the bulk. 

A final scenario for the discrepancy is that the band 
structure is distorted near the surface due to space-charge 
accumulation. This is known to occur in many semicon- 
ductors, such as InSb or CdTe 15, whereby the bulk 
band structure bends as the surface is approached. Typ- 
ically, such bending occurs over a surface depletion layer 
Zd, which can be calculated by solving the Poisson equa- 



tion to yield 



wnAy/en. [171] ■ where k is the DC 



dielectr ic p ermittivity (estimated from these samples as 
~ 113 [18])and Ay is the difference in energy between 
the surface and bulk state. We estimate z^ 1 ~ 40nm for 
the low carrier density samples and zf 4 ~ 2nm for the 
high carrier density samples. A schematic representation 
of the band bending is shown in Figure [5] The present 
argument suggests that discrepancies between ARPES 
and SdH can be explained, even expected for low carrier 
density samples. In addition, due to the small value of 
n, these samples are likely more susceptible to a small 
amounts of surface contamination, especially if the un- 
contaminated surface Ep is in the gap (as illustrated in 
by ARPES on atmosphere exposed samples). This may 
help explain why there is some variability in the Fermi 
level of ARPES data but not in the SdH data. 

Much theoretical work has emerged on the dra- 
matic con seq uences of the surface state on transport 
properties [19j. Yet over several decades of experimental 



study, such properties have not been observed. Recently, 
Aharanov-Bohm and universal conductance fluctuations 
have been observed which may be due to the surface state 
nj ElD], but even in these cases the conductance appears 
bulk at the temperatures considered. Conventionally, 
such intrinsically doped materials can become 'insulators' 
by either a Mott-like or an Anderson transition. The first 
can occur when the Bohr radius = nh 2 /m*e 2 falls 
below the Thomas- Fermi screening length Xtf, so that 
wavefunctions cannot overlap. This can be estimated us- 
ing A^ F = ne^/ (2ire 2 g(Ep)), where g{Ep) is the number 
of states per unit volume per unit energy, estimated by 
Middendorff et aZ.[2(|. In the present case, the large k 
and small m* tend to make as very large. For the low- 
est carrier density samples investigated here ~ 3nm 
and Xtf ~ 4nm, which places is this material on the 
metal-insulator boundary. The carrier density can also 
be reduced by introducing foreign dopants which 'drain' 
the excess carriers and pin the chemical potential /i in 
the gap. For hydrogenic like impurities this can be very 
effective, but in the present materials impurity bands of- 
ten form instead. At high enough impurity densities the 
carriers may become Anderson localized. Such samples 
are characterized by a high carrier density with very low 
mobility, leading to a negative gradient in the tempera- 
ture dependence of the resistivity. This may be the case 
for example in BLjShyP^Ses which has p ~ 30mf2cm 



2l| . An Anderson 



yet a carrier density n e ~5xl0 18 cm~ 3 
insulator is generally bad news for topological insulators, 
because even though at zero temperature the bulk con- 
ductivity (7 = 0, at finite temperature the transport may 
remain dominated by bulk hopping mechanisms. 

In conclusion, the present study reveals substantial 
agreement between transport and ARPES measurements 
of the Fermiology of Bi2Se3, in particular for samples 
with large carrier densities. However, for samples with 
carrier densities approaching 10 17 cm~ 3 , discrepancies 
emerge as to the exact position of the Fermi level. We 
have confirmed the bulk nature of the transport by the 
thickness dependence of the Hall effect, resistivity and 
mobility. Furthermore SdH data is highly consistent be- 
tween different samples from the same batch. Interest- 
ingly, the carrier densities measured here are an order 
of magnitude smaller than those of the topolo gica l in- 



sulators recently reported in the literature l^lSt M, EH- 
ARPES and STM have been invaluable tools in reveal- 
ing the physics of topological insulators, providing com- 
pelling evidence for the presence of the topologically pro- 
tected Dirac surface state. The present results should 
stimulate further theoretical work as to the consequences 
of the coexistence of bulk and surface states in a single 
sample as well as innovation in novel ways to fabricate 
these materials so the bulk state can be cleanly elimi- 
nated. 
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